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ABSTRACT 

Context. The common flare scenario comprises an acceleration site in the corona and particle transport to the chromosphere. Using 
satellites available to date it has become possible to distinguish between the two processes of acceleration and transport, and study 
the particle propagation in flare loops in detail, as well as complete comparisons with theoretical predictions. 

Aims. We complete a quantitative comparison between flare hard X-ray spectra observed by RHESSI and theoretical predictions. This 
enables acceleration to be distinguished from transport and the nature of transport effects to be explored. 

Methods. Data acquired by the RHESSI satellite were analyzed using full sun spectroscopy as well as imaging spectroscopy methods. 
Coronal source and footpoint spectra of well observed limb events were analyzed and quantitatively compared to theoretical predic- 
tions. New concepts are introduced to existing models to resolve discrepancies between observations and predictions. 
Results. The standard thin-thick target solar flare model cannot explain the observations of all events. In the events presented here, 
propagation effects in the form of non-collisional energy loss are of importance to explain the observations. We demonstrate that those 
energy losses can be interpreted in terms of an electric field in the flare loop. One event seems consistent with particle propagation or 
acceleration in lower than average density in the coronal source. 

Conclusions. We find observational evidence for an electric field in flare loops caused by return currents. 
Key words. Sun: flares - Sun: X-rays, gamma-rays - Acceleration of particles 



1. Introduction 

Solar flares have been studied in detail both observation a lly an d 
theoretically ever since their discovery by Carrin gton (1859). 
Although increasingly more sophisticated instrumentation pro- 
vides ever more detailed data, we still lack the basic understand- 
ing of many processes at work in a solar flare. 

The common flare picture as deduced from hard X-ray 
(HXR) observations feat ures an HXR source in t he corona (coro - 
nal or loop-top source, iFrost & Dennii 119711: iHudsonl Il978h . 
and tw o or more HXR so urces (footpoints) in the chromo- 
sphere (Hov nget alj 1 1 98 lb . These sources are thought to be 
due to bremsstrahlung emission produced by fast electrons ac- 
celerated somewhere above the loop. If we assume that a sin- 
gle particle beam creates both coronal and footpoint emis- 
sion, the most basic model would involve thin target emis- 
sion at the top of the coronal loop and thick target emission 
from the footpoints, which b oth produce characteristic spectra. 
IWheatland & Melrose! JT995) developed a more sophisticated 
model (intermediate thin-thick target or ITTT model) to fit ob- 
ser vations by Yohkoh. T hey based their model on observations 
by iFeldman et alJ {1994), who found high column densities at 
the loop top, which might act as a thick target below a certain 
electron energy. In the ITTT model, the shape of the coronal and 
footpoint non-thermal spectra and the relation between them, ob- 
served by Yohkoh, can be explained. The column density in the 
coronal source determines a critical energy level for the elec- 
trons. Electrons that have an energy below this critical energy 
are stopped in the coronal region. Consequently, the distribu- 
tion of electron energies measured at the footpoints is depleted 
in low energy electrons. If the column density is high, the coro- 
nal source may act as a thick target to electrons of energies as 
high as 60 keV, which would leave almost no footpoint emis- 



sion. Observational evidence for su ch coronal thick targets were 
found in RHESSI observations (eg. lVeronig & Brownll2004l) . 

Less extreme cases, flares with one or more footpoints, have 
frequently been observed by RHESSI. To study the spectral time 
evolution of individual sources , five w ell-observed events were 
analyzed by Battaglia & Bend (H006) who focused on the dif- 
ferences between the spectral indices of coronal and footpoint 
spectra. They found that in two of those events, the differences 
at specific times as well as the time-averaged difference was sig- 
nificantly larger than two, ruling out a simple thin-thick target 
interpretation. In lBattaglia & Bena (120071) . the spectra of the five 
events were compared with the predictions of the ITTT model. 
The authors exploited the order of magnitude improvement in 
spectral resolution of RHESSI over the 4-point Yohkoh spectra 
and showed that most RHESSI observations could not be ex- 
plained by the ITTT model. 



Battaglia & Bend (120071) proposed that by considering non- 
collisional energy loss inside the loop this inconsistency could 
be resolved. A possible mechanism that causes non-collisional 
energy loss is an electric field. Accelerating electrons out of 
the coronal source region drives a return current to main- 
tain charge neutrality in the whole loop. For finite con- 
ductivity, Ohm's law implies that an electric field must be 
present. The beam electrons lose energy because of work ex- 
pended in moving inside the electric potential. This produces 
a change in the shape of the electron spectrum at the foot- 
points. The formation and evolution of these return currents 
were studied by various authors (e.g. Knight & Sturrockl 1977 : 
ISpicer & Sudani Il984t Larosa & Emslid Il989t Ivan den Oordl 
1990). IZharkova & Gordovsk yy (2006) proposed that return cur- 
rents could explain the high energy break observed in flare HXR 
spectra. Most studies have, however, been theoretical proposals 
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or numerical simulations, based on standard flare values, that do 
not at tempt to explain or rep roduce true solar flare observations. 

IBattaglia & Benzl d2007l) compared RHES SI spectra to the 
ITTT-model of IWheatland & Melrosel d 1995b . demonstrating 
that the qualitative shape and relations between coronal and 
footpoint-spectrum often do not agree with the model predic- 
tions. In this study, we take an additional step by completing a 
quantitative analysis of the relation linking coronal and footpoint 
spectra in the context of the thin-thick target model; we demon- 
strate that, in some cases, electric fields related to return currents 
can indeed explain the relation between coronal and footpoint 
spectra. 

In Sect. [2] we summarize the basic physical concepts applied 
in the paper. Section[3]provides a brief overview of the analyzed 
events and a description of the spectral analysis. In Sect. [4] we 
describe our calculation of the energy loss required to repro- 
duce the observed footpoint spectrum, constrained by the coro- 
nal emission. Our results are presented in Sect. [5] In Sect. [6] we 
link those results to the concept of return currents. 

2. Thin and Thick target emission 

Two types of bremsstrahlung emission are distinguished. If the 
electrons pass a target without losing a significant amount of 
ener gy, the correspondin g emission is referred to as thin tar- 
get dDaflowe & Lin|[l973l) . This situation is expected to occur 
in coronal regions when electrons pass through a target of in- 
sufficient column density to stop them. If the electrons are fully 
stopped inside t he target, the resultant emission is called thick 
target emission dBrown|[l97ll) . This is the case for the dense 
chromospheric material at the footpoints. 

For an input power-law electron distribution of the shape 
F(E) = AeE~ s , the non-relativistic bremsstrahlung theory pre- 
dicts power-law photon spectra 



1(e) oc e 1 where 



y = S + 1 in the thin target case 
y = S - 1 in the thick target case 



The observable distinction between the two emission mecha- 
nisms is a difference Ay of value 2 in their observed photon 
spectral indices. 

Assuming a final column density in the coronal source, the 
coronal source spectrum is a thick target at low energies and 
a thin target at high energies with a break at some critical en- 
ergy. The footpoint spectrum is depleted at low energies, as 
low energy electrons do not r each the chromosphere . An ill us- 
tration of this can be fo und in Wheatl and & Melrosel d 1995b or 
iBattaglia & Benzl d2007l) . 

In the events that we analyze here, a thermal component is 
present in all observations. Observation of the non-thermal emis- 
sion is therefore only possible at photon energies higher than 15 
keV. For these energies we can assume that the coronal source is 
a pure thin target and the footpoints are a pure thick target. 

3. Event description and spectral analysis 

The two events that we analy ze were described in detail by 
IBattaglia & Benzl d2006l 120071) . They were selected because 
of the significant differences between the footpoint and coro- 
nal non-thermal spectral index. The first event occurred on 24 
October 2003 around 02:00 UT (GOES M7.7), the second on 
13 July 2005 around 14:15 UT (GOES M5.1). Both events oc- 
curred close to the limb but were not occulted; two footpoints 
were therefore fully observed in both cases. RHESSI light curves 



from the time of main emission are shown in Fig.Q] An EIT im- 
age of the 24 October 2003 event and a GOES SXI image for the 
event of 13 July 2005 are presented for orientation. The contours 
of the coronal source and the footpoints from RHESSI images 
are overlaid. 



3. 1 . Spectral fitting and analysis 

The two events were a nalyzed using imaging spectroscopy wit h 
the PIXON algorithm dMetcalf et al.lll996tlHurford et alj|2002l) . 
Images were made in a 30 second time interval during which the 
flux was sufficiently high for good images but pile-up was low. 
The image times are given in Tableland shown in Fig.Q] The 
spectra of the footpoints and the coronal source were measured 
and fitted. The regions of interest from which the spectra were 
computed were chosen to be a circle around the coronal source 
and a polygon around the footpoints to include all of the emis- 
sion at all energi es. The effects of this m ethod of region selection 
are discussed in Battagli a" et al.l d2005l) . As a simplification, the 
footpoints were treated as one region and the spectrum was fit- 
ted with a single power-law. In the presence of the footpoints, 
the non-thermal emission in coronal sources is difficult to ob- 
serve. Therefore, two methods of fitting the coronal source were 
used. First, a thermal component was fitted to the spectrum at 
low energies and a single power-law to the energies higher than 
about 25 keV. As a second method, the ful l sun t hermal spec- 
trum was fitted. As shown in Battaglia et alJ d2005l) . the thermal 
emission observed in full sun spectra is mostly coronal emission. 
We therefore used the thermal full sun fit as an approximation 
to the coronal thermal component and completed a power-law 
fit at the higher energies, while the thermal emission was fixed. 
This supports the idea that non-thermal emission exists in the 
coronal source and provides an estimate of the accuracy of the 
non-thermal coronal fit. The energy ranges for the fits were 8-36 
keV for the coronal source and 24-80 keV for the footpoints. All 
fit parameters are provided in Table Q] 

In the thin-thick target model an electron beam is assumed 
to be injected into the center of the coronal source. The column 
depth that the electrons travel through in the corona is then 
AN = n e ■ I, where the path length I is half the coronal source 
length. From RHESSI images in the 10-12 keV band, the source 
area A was measured to be the 50% contour of the maximum 
emission. We approximate the source volume to be V = A 3 ^ 2 
and the path length to be / = yA/2. Using the observed 
emission measure EM, the particle density is computed to 
be n e — y/EM/V which corresponds to a column depth of 
AN = y/EMA~ l l 4 /2 expressed in observable terms. A volume 
filling factor of 1 was assumed for the computation of the 
density, which will be improved in Sect. 15.21 The emission 
measures were taken from the spectral fits to the coronal source 
and to full sun spectra. Additionally, temperatures and emission 
measures observed by GOES were included. This provides a 
range for the emission measures, temperatures, and column 
depths, and an estimate of their uncertainty. 



4. Method 

Starting from the assumption that the observed coronal spec- 
trum at high photon energies is caused by thin target emission, 
we compute the electron distribution and therefore the expected 
footpoint photon spectrum. This is completed via the following 
steps. 
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24 October 2003 13 July 2005 




02:20 02:30 02:40 02:50 03:00 03:10 03:20 14:10 14:15 14:20 14:25 14:30 14:35 14:40 
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Fig. 1. Top: RHESSI light curves in the 6-12 and 25-50 keV energy band. The analyzed time interval is indicated by the vertical 
bar. Bottom: SOHO/EIT image at 24 October 2003 02:47:31 {left), GOES SXI image at 13 July 2005 14:19:04 (right). The 30 %, 
50 % and 70 % contours from RHESSI Pixon images of the coronal source at 12-16 keV {solid contours) and 20-24 keV (dashed 
contours, mainly non-thermal emission) are overlaid along with the 50% contour of the footpoint sources at 25-50 keV. 



1 . We assume that the observed coronal photon spectrum can 
be fitted by a power law: 

F c : bs ie)=K s ^\ (1) 

where A" is the normalization and y cs the photon spectral 
index. 

2. Using thin target emission, the injected electron spectrum 
F(E) is then proportional to 

F(E) = A E E- 5 ~ %E-v" +1 (2) 

where AN is the column depth the electrons t ravel through 
inside the coronal target (Dat lowe & Lin|[l973l) . 

3. The expected thick target emission F^ p (e) caused by this 
electron distribution in the footpoints can be computed as 
follows dBrownlll97ll) : 

r A cs 

Fl P xp (e)=At Kp e-^>~^e-*"-V (3) 

The superscripts fp and cs denote the footpoint and coronal 
source values, respectively. 

4. The normalization and spectral index of F~£p{e) is compared 
to the observed footpoint spectrum F^f ' (e). 

In thin-thick target models, the difference in spectral index Ay = 
\y cs - y' p \ is 2. As the observed difference is larger than 2 in 
the selected events, a mechanism has to be found that causes 



the electron spectrum to harden while the beam passes down 
the loop. We present a mechanism that assumes an electric field 
which causes electrons to lose the energy £/ oss independently of 
the initial electron energy. The resulting spectrum is flatter, al- 
though not strictly a power-law function anymore (Fig. [2^). The 
deviation becomes substantial below 2 &i oss . 

The necessary energy loss is determined as follows: 

1 . We start with the coronal electron distribution as found from 
Point|2]in the above list. 

2. By assuming a thin target, the electron distribution leaves 
the coronal source and propagates down the loop. A constant 
energy loss &i oss is subtracted from the electron energies as 
the energy loss is independent of the electron energy. 

3. We compute the expected thick target photon spectrum 
Fexp( £ ) from this altered electron spectrum. 

4. A power law is fitted to F^ p (e). The fitted energy range is 
30-80 keV. This is the range for which footpoint emission is 
typically observed. 

The relation between the energy loss experienced and the corre- 
sponding photon spectral index of the best-fit power law func- 
tion depends on the initial electron spectral index 5 and the en- 
ergy loss Si oss . It is equivalent to the elementary charge times 
the electric potential between the coronal source and the foot- 
points. If the initial electron spectral index is 8 for instance, the 
thick target photon spectral index without energy loss is 7. With 
increasing energy loss, this value decreases rapidly. The effect 
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is less pronounced when the initial electron spectrum is harder. 
This is shown in Fig.[2j3 for several values of 6 and £/ oss . Using 
the curves in this figure, we can easily determine the energy loss 
that causes an electron spectrum of spectral index 5 to result in a 
fitted photon spectral index yf p . 

5. Results 

5.1. Observed spectra 

Figure [3] shows the observed spectra overlaid with the spectral 
fits. As indicated in Sect. 13. II the thermal fits differ slightly from 
each-other. The main reason why the fits do not agree is the 
wider energy binning adopted by imaging spectroscopy. With 
this binning, the atomic lines are not resolved, contrary to full 
sun spectroscopy. 

Using the different fitting methods as an estimate of the 
uncertainty, an average difference in spectral index of Ay = 
3.55±0.07forthe event of 24 October 2003 and Ay = 2.45+0.35 
for the event of 1 3 July 2005 is found between the coronal source 
and footpoints. 

5.2. Expected footpoint emission and energy loss 

As described in Sect. @] we computed the electron distribution 
from the coronal source photon spectrum and the expected thick 
target emission (footpoint spectrum) caused by this electron dis- 
tribution. Figure @] shows the measured spectra, the expected 
footpoint spectrum from a pure thick target, and the footpoint 
spectrum when introducing energy loss. 

To compute the electron flux, the coronal column density is 
required. As given in Table Q] three different values of the col- 
umn density were estimated. Using those, we are able to repro- 
duce a range of possible electron spectra (Tabled and, therefore 
footpoint spectra. The confidence range of the footpoint spectra 
is indicated by a light gray (green) area in Fig. [4] 

The derived energy loss depends on the fitted coronal and 
footpoint spectra. For the two different coronal fitting methods, 
a range of loss energies &i oss = [58.0,59.4] keV is found for 
the event of 24 October 2003 and S hss = [8.7, 26] keV for the 
event of 13 July 2005. The normalization of the new spectrum 
depends on the initial electron distribution. The initial electron 
distribution is computed according to Eq. ©. It depends on the 
column depth. If the column depth is lower, more electrons are 
needed to produce the same X-ray intensity. The resulting range 
of possible spectra is shaded in dark-gray in Fig. [4] As shown 
in the figure, the footpoint spectrum with energy loss reproduces 
well the observed footpoint spectrum for the event of 13 July 
2005. 

During the event of 24 October 2003, the predicted footpoint 
spectrum is, however, an order of magnitude less intense than 
observed. In the context of the ITTT model, this implies that the 
electron flux density emanating from the coronal region is higher 
than predicted. The discrepancy may be explained by density in- 
homogeneities in the coronal source resulting in a smaller effec- 
tive column density. The electron flux is underestimated if the 
non-thermal X-ray emission originates in regions that are less 
dense than average. In the following we therefore assume that 
the coronal source has an inhomogeneous density; this is rep- 
resented by dense regions with filling factor smaller than 1 for 
thermal emission and, for the non-thermal coronal source in the 
24 October 2003 event, a density that is lower by an order of 
magnitude. The observations do not allow to determine the fill- 
ing factor. As can be seen from Fig.Q] the coronal source in the 



event of 13 July 2005 is very compact, while the source in the 
event of 24 October 2003 is more extended, showing isolated 
intense regions. This supports the assumption that the density 
is inhomogeneous in the 24 October 2003 coronal source and 
the true column density in the X-ray emitting plasma might be 
smaller than deduced from the measurements. Since F^ p (e) is 
proportional to 1 / AN, an effective column density of an order of 
magnitude less than the observed could produce the observed 
footpoint spectrum. In the computations presented in Sect. [6] 
we assume an effective column density AN e ff = AAV 14 f° r the 
event of 24 October 2003. 



6. Return current and electric field 

In the above analysis, we assumed that the electrons experienced 
a constant energy loss while streaming down the loop. We now 
demonstrate that this energy loss could be caused by an elec- 
tric potential in the loop that drives a return current. There was 
much controversy surrounding the pre cise physical mechanism 
that generates the return current (eg. iKnight & Sturrock1ll977l: 
Spicer & Sudani [19841: Ivan den Oordlll990l) . The basic scenario 
is the following: We assume that the electrons are accelerated 
in the coronal source region. When a beam of accelerated elec- 
trons, which is not balanced by an equal beam of ions, leaves 
this region, a return current prohibits charge build-up and the 
induction of a beam-associated magnetic field. In the return cur- 
rent, thermal electrons move towards the coronal source. Since 
their velocity is relatively small, they collide with background 
ions and cause resistivity. Ohm's law then implies the presence 
of an electric field in the downward direction. The return current 
density j ret can be de rived from t he equation of motion for the 
background electrons (Benz 2 0021) . 



dv e e 

— + (v • V)v = E ind (v X B) - v e>i v 

ot m mc 



(4) 



where Ei„d is the electric field induced by the return current, v 
is the mean velocity of the electrons that represent the return 
current, and v e .,■ is the electron-ion collision frequency. Using 
j re t - -env, this expression can be written as 



d 



(wp 2 e 

jret = — ; Ei n d H \jret X (^0 + Bj„d)], 

4tt mc 



(5) 



where Bo is the guiding magnetic field and B m( i the field induced 
by the beam. We neglect the last term on the right side of Eq. 
by assuming that the beam and return currents are anti-parallel, 
oriented along the guiding magnetic field, and that the perpen- 
dicular component of B,„</ vanishes. We then obtain 



d 

— + y e 

dt 



Jret — 



An 



J ind ■ 



(6) 



Since we consider a fixed time interval that is far longer than 
the collision time, we assume a steady state, neglecting the time 
derivative of the return current. The equation then takes the form 
of the classical Ohm's law: 



Jret — 



Ei„d = o~E, 



ind- 



(7) 



We estimate whether the energy loss computed in Sect. [5] 
is caused by this electric field. From the observations, we es- 
timated the electron loss energy &i oss that the electrons expe- 
rience in the loop (Sect. |5.2| >. Assuming this loss is caused by 
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Fig. 2. a) Normalized electron power-law spectrum (solid) and altered spectrum due to a constant energy loss of 30 keV (dashed). 
b) Relation between loss energy fi/ OSJ and fitted thick target photon power-law spectral index yf p for initial electron spectral index 
6 = 3, 4, 5, 6, 7, 8 in the accelerator. 



Table 1. Overview of main event properties and fit parameters 



Time interval 


24 October 2003 02:48:20-02:48:50 


13 July 2005 14:15:00-14:15:30 


Area [cm 2 ] / Volume [cm 3 ] 


7.9 


10 18 /2.2 


10 JS 


1.7 • 


10 18 /2.2- 


10 2/ 




full sun 


imspec 


GOES 


full sun 


imspec 


GOES 


Temperature [MK|^ 


21.6 


23.2 


15.4 


23.8 


22.3 


18.1 


Emission measure [10 49 cnT 3 ] 


0.98 


0.44 


3.4 


0.22 


0.22 


0.47 


Electron density [10 10 cnT 3 ] 


2.1 


1.4 


3.9 


3.2 


3.2 


4.6 


Column density [10 19 cirr 2 ] 


2.9 


2.0 


5.5 


2.1 


2.1 


3.0 




footpoints 


cs fit 1 


cs fit2 


footpoints 


cs fit 1 


cs fit 2 


y- 


2.6 


6.2 


6.1 


2.9 


5.1 


5.6 


F50 [photons cm 2 s 'keV '] 


2.0 


0.07 


0.07 


0.88 


0.06 


0.05 



" Thermal parameters for three different measuring methods (RHESSI full sun fit, RHESSI imaging spectroscopy, GOES). See also comment 
in Sectl5~T1 

4 Two different values (cs fit 1, cs fit 2) for the non-thermal coronal fit, distinguishing the two fitting methods used (compare Sect |3.1t . 




10 100 10 100 

Energy [keV] Energy [keV] 



Fig. 3. Coronal and footpoint source spectra overlaid with the according fits. Dots are the measured footpoint spectrum, the dashed- 
dotted line indicates the fit to this spectrum in the range 30-80 keV. Squares indicate the observed coronal source spectrum. The 
solid lines provide the thermal and non-thermal fits as found from imaging spectroscopy. The dotted lines give the thermal fit to the 
full sun spectrum and the resulting non-thermal fit in imaging spectroscopy. 
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Fig. 4. Observed (fitted) power-laws of the non-thermal coronal source and footpoints (solid). The light-gray (green) area indicates 
the range of expected footpoint spectra without energy loss. The dark-gray (red) area marks the range of expected footpoint spectra 
when energy loss is applied to the electrons to find the same spectral index as the observed footpoint spectrum. The dash-dotted lines 
give the expected spectra at the footpoints if the only transport effect was Coulomb collisions of the beam electrons (cf. Sect. I7.4I >. 



the induced electric field £,„,/ and across the distance s from the 
coronal source to the footpoints (i.e. half the loop length), we 
compute the electric field to be 

Em = — —■ (8) 
e ■ s 

Using Spitzer conductivity (Spitzerl [l965l) . the term for the 
return current is related to the observed loss energy by: 

jre, = 6.9 ■ 10 6 r^^ [statamp/cm 2 ], (9) 

where Ti oop is the temperature in the loop. 

On the other hand, the beam current density can be written 

as 

jbeam = F ""^ ■ e [statamp/cm 2 ], (10) 

A fp 

where A f p is the total footpoint area. The total electron flux per 
second F'°'(E) is computed from the observed electron spectrum 
as follows: Let the electron spectrum be F(E) - A e E~ s . The total 
flux of streaming electrons per second above a cutoff energy E cut 
is then: 

F m (E)= F(E)dE=—^-E^- 1) . (11) 
Je„„ o - 1 

In a steady state, the relation 

Jbeam ~ J ret 

(12) 

is valid. 

Comparing the beam current as described in Eq. (fTOb with 
the return current from the observed energy loss according to 
Eq. (O, we test whether the assumption of Spitzer conductivity 
holds. 



6.1. Results 

Table [2] presents the relevant physical parameters necessary for 
the derivation of the beam- and return currents. For Spitzer con- 
ductivity, the loop temperature Ti oop is required. Its is expected 
to have a value between the coronal source temperature and the 
footpoint temperature (see Table |2|. As a first assumption, a 
mean temperature of Ti oop = 15 MK is chosen. The loop length 
is evaluated from RHESSI images, approximating the distance 
between the sources from the centroid positions and assuming 
a symmetrical loop structure. This provides a typical half loop 
length of 4 • 10 9 cm. The footpoint area is measured from the 
50% contour in RHESSI images in the 25-50 keV energy range, 
yielding a total footpoint area of « (6 - 7) • 10 17 cm 2 . The beam 
current density depends critically on the electron cut off energy 
E cut . We use a value of 20 keV. This gives an approximate lower 
limit to the total amount of streaming electrons. 

Using the presented observations, Eq.dTOl) and by assuming 
Spitzer conductivity (Eq.[9]), the return current results to be of an 
order of magnitude higher than the beam current. This contra- 
dicts the assumptions of a steady state, and is also unphysical. 



7. Discussion 

7.1. Instability 

In Sect. I6.ll we assumed Spitzer conductivity when comput- 
ing the return current which produced the unphysical result of 
jrei > jbeam- Using Eq. ®, the loop temperature required to 
maintain equality between the return current and beam current 
Uret = jbeam) can be computed. In the 24 October 2003 event, the 
loop temperature Ti oop would need to be smaller than 10 MK; in 
the 13 July 2005 event, Ti oop should be less than 3.9 MK. Such 
low loop temperatures are highly unlikely. However, it is pos- 
sible that the return current is unstable to wave growth. For an 
extende d discussion of instabilities in parallel electric currents, 
see e.g. [Benz (2002). Instability causes an enhanced effective 
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Event 



24 October 2003 



13 July 2005 



Assumed loop temperature T [oop [MK] 
1/2 Loop length s [cm] 
Electron flux F(E) [s~' keV _1 ] 
Electron cutoff energy [keV] 
Total footpoint area [cm 2 ] 
Eioss [keV] 

Electric field strength [statvolt/cm] 
[statamp I cm 1 ] 



6.33-10 41 £ 



15 

3.2 • 10 9 
"-1.5-10rf 42 £- 52 
20 

7.2 • 10 17 
58-59.4 
(6 - 6.3) • 10- 8 



•2.5) • lfF" 
14.4) • 10 9 



15 

4.3 • 10" 

1.7-10 40 £- 41 - 8.4- 10 40 £- 4 - 6 
20 

6.2- 10 17 
8.7-26 
(6.7 - 20.3) • 10" 9 



10 8 



J beam 



[statamp /cm-] 



(2.4 
(5.1- 



(2.7- 
(1.9- 



8.1) 
3.6) 



collision frequency of electrons in the return current and there- 
fore a lower effective conductivity. The ion cyclotron instability 
develops if the drift velocity of the beam particles Vd exceeds 



the thermal ion velocity v'™ as follows: 



(13) 



with v"!" = A /^ s - J and T e and T: being the electron and ion tem- 

th y mi L ' & 

peratures, respectively. 

We assume a steady state for which jbeam = jret = n e eVd, 
where Vd is the mean drift velocity of the electrons constituting 
the return current. We therefore express Vd as 



V d 



Jbe, 



(14) 



and substitute this expression and that for vj?" in Eq. ( TT3T >. 
Assuming T e = Tj = Ti oop and solving Eq. (fl3b for Ti oop the 
instability condition holds 



loop 



< 2.3 ■ 10" 



8 / Jbe. 



[K]. 



(15) 



Since the loop temperature and density are not known ex- 
actly, this relation is illustrated in Fig.|5]for several values of n e 
and Tioop typical in flare loops. For the values of jbeam found for 
the observations of the two events, we find that instability oc- 
curs in the 24 October 2003 event for all values of n e and Ti oop 
in Fig. [5] For the 13 July 2005 flare, three distinct regions in the 
diagram can be found. The solid line indicates the relation of 
Eq. ( fT3T ). Below this line, the return current is unstable. At high 
densities and low temperatures (lower right), the return current 
is stable and jbeam = jret with Spitzer conductivity. The range of 
beam currents jbeam deduced from the data allows for loop tem- 
peratures < 3.9 MK. In the upper right quadrangle, Spitzer con- 
ductivity would imply j ret > jbeam, which is unphysical. If the 
loop was in this parameter range, the current instability would 
be most likely saturated and T e > Tj. This would shift the insta- 
bility threshold in Fig. [5] to the right. Further, a loop in the state 
presented by the uppermost part of the figure (temperature above 
10 MK, high density) would be detectable by the RHESSI satel- 
lite even in the presence of the coronal source. Since no loop 
emission is observed, we conclude that the loop is either less 
dense, cooler or both. The values in the upper right quadrangle 
are therefore unlikely. 
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Fig. 5. Region of instability in the density/temperature space for 
the event of 13 July 2005. The grey region indicates the densities 
and temperatures Ti oop for which the return current is unstable 
(Eq.[T3ll. In the lower right part, the return current is stable and 
jbeam = jret with Spitzer conductivity. The current in the event 
of 24 October 2003 is unstable for all values of density and tem- 
perature in the Figure. 



which the thermal and non-thermal components of the spectrum 
intersect. Values around 20 keV or higher are also supported by 
detail ed studies of the exact dete r mination of low-energy cutoffs 
(eg. ISaint-Hilaire & Benzl 120051: IVeronig et alj l2005t [Sui et alj 
2007). What if the cutoff energy were substantially lower than 
20 keV? A cutoff energy as low as 10 keV would increase the 
total electron flux and therefore the beam current by an order 
of magnitude, leading to j re , ~ jbeam for Spitzer conductivity. 
Conductivity could then not be reduced significantly by wave 
turbulence, and instability would be marginal. If the low en- 
ergy cutoff were even lower than 10 keV, we would find that 
jret < jbeam- This could not be explained in terms of the model 
used here. 



7.2. Low energy electron cutoff 

In the above computations, a value of 20 ke V for the electron cut- 
off energy E cut was assumed. This value is within the range for 



7.3. Source inhomogeneity and filling factor 

In the above paragraphs, it was demonstrated that the energy 
loss for the electrons due to an electric field could resolve the 
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inconsistency in the difference between footpoint and coronal 
source spectral indices (Sect. |5]l- F° r the event of the 24 October 
2003, this produces footpoint emission that is lower than the ob- 
served emission (Fig.|4]i. Images show that the coronal source in 
this event is not compact, but extended with brighter and darker 
regions. It is therefore possible that the standard density esti- 
mate, which favors high density regions, produces higher den- 
sities than average and that the effective column density of the 
regions, where the largest part of the non-thermal emission orig- 
inates, is lower. This would provide a higher expected footpoint 
emission, in closer agreement with observations. 

As mentioned in Sect. 13.11 a filling factor of 1 was used for 
the computation of the column density. A filling factor smaller 
than 1 would lead to even higher densities of the SXR emitting 
plasma. However, this would not affect the lower effective col- 
umn density of the regions, where the HXR emission originates 
when assuming source inhomogeneity. 

7.4. Collisions and other possible scenarios 

While return currents may not be the only means of attaining 
non-collisional energy loss, they are the most obvious and best 
studied. However, other scenarios are conceivable, which could 
produce a harder footpoint spectrum (or a softer coronal spec- 
trum). In the model presented here, collisional energy loss of 
beam electrons is neglected. This is a valid assumption for the 
following reasons: If collisions of the beam electrons in the 
loop were to play an important role, significant HXR emission 
should originate in the loop. Within the dynamical range limita- 
tions of RHESSI, this is not the case. However, GOES SXI and 
SOHO/EIT images imply that the loop is filled with hot material. 
To study possible effects of collisions, we compared the change 
in the electron spectrum and the resulting footpoint spectrum for 
collisional energy loss and energy loss due to the electric field. 
The change in the electron spectrum due to collisions depends on 
the colum n depth through wh i ch th e beam passes and was com- 
puted by iLeach & Petrosianl (Il98ll) and iBrown & McClvmontl 
(Il975h . We assume a column depth derived from the density in 
the loop times half the loop length. Assuming the same density 
as in the coronal source, we derive an upper limit to the col- 
lisional effects. The expected footpoint spectrum from purely 
collisional losses is indicated in Fig. [4] as dash-dotted lines. 
Collisions affect the low energetic electrons most where a sig- 
nificant change in the spectrum is found. At the higher ener- 
gies observed in this study, the spectrum does not change sig- 
nificantly. The neglect o f collisional effects is therefore justified. 
IBrown & Mallikl d2008) showed that in certain cases, emission 
from non-thermal recombination can be important, generating a 
coronal spectrum that is steeper than expected by the thin-target 
model. This could also produce a difference in the spectral index 
that is larger than two. Acce leration over an extended region (as 
proposed bv lXu et aL 2008) could alter the electron distribution 
at the footpoints. If the distribution was harder at the edge of the 
region, a spectral index difference larger than 2 would result. A 
thorough comparison of such models with observations may be 
the scope of future work. 

8. Conclusions 

The spectral relations between coronal and footpoint HXR- 
sources provide information about electron transport processes 
in the coronal loop between the coronal source and the foot- 
points. Most models neglect these processes in the prediction 
of the shape and quantitative differences between the source 



spectra. As shown by Battaglia & Benz (120071) . the observations 
of some solar flares do not fit the predictions of such mod- 
els, in particular the i nterm ediate thin-thick target model by 
IWheatland & Melrose! (TT995): there is a discrepancy concerning 
the difference in coronal and footpoint spectral indices, which is 
expected to be 2. 

We have analyzed the two out of five events that display a 
spectral index difference larger than two in more detail. Such a 
behavior can be attributed to energy loss during transport that is 
not proportional to electron energy, but Si oss /E is larger at low 
energies. Such an energy loss causes the footpoint spectrum to 
flatten, which increases the difference in spectral indices. Two 
loss mechanisms come to mind immediately: Coulomb colli- 
sions and an electric potential. Figure 4 demonstrates that the 
assumption of an electric potential reproduces the observations 
more accurately. 

In one of the two events, there remains a discrepancy be- 
tween the observed and expected footpoint emission, such that 
the electron flux at the footpoints is larger than predicted. This 
flux was estimated from the observed non-thermal HXR (pho- 
ton) flux and the observed thermal emission of the coronal 
source. We attribute the discrepancy to propagation or accelera- 
tion in low density plasma, which also heats the adjacent high- 
density regions. 

The energy loss can therefore be explained by an electric 
field in the loop associated to the return current, which builds up 
as a reaction to the electrons streaming down the loop and the as- 
sociated beam current. In a steady state (jbeam = jret), the return 
current is unstable to wave growth in one event for all realis- 
tic temperature and density parameters in the loop. The kinetic 
current instability drives a wave turbulence that enhances the 
electric resistivity by many orders of magnitude. This anoma- 
lous resistivity in turn significantly enhances the electric field. 
In the event of 13 July 2005, the return current may be stable if 
the loop density is high and the temperature is low, and Spitzer 
conductivity is applied. Both cases (out of five) present strong 
evidence for a return current in flares for the first time. 

Transport effects by return currents constitute a considerable 
energy input by Ohmic heating into the loop outside the accel- 
eration region. It may be observable in EUV. Comprehensive 
MHD modeling including the coronal source, the footpoints, 
and the region in-between, may be the goal of future theoreti- 
cal work. 
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